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ABSTBACT 



Wave spectra calculated using the Parameterized 
Nonlinear Wave Solution developed by Le Mehaute et al. 
(1984) are compared with field data acquired at Leadbetter 
beach , Santa Barbara, California. The parameterized solu- 
tion satisfies the nonlinear free surface boundary condi- 
tions to a specified degree of accuracy and is expressed in 
terms of a converging truncated Fourier series. The wave- 
number, surface profile and wave orbital velocities are 
determined by the wave height and wave period at the local 
depth of water. Spectral components are compared between 
the model results and field data. Good agreement is observed 
for waves corresponding to Drsell numbers ranging from 25 to 
75. For large Ursell numbers (strong nonlinear effects) the 
parameterized model underestimates the data. 
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I. INTRODUCTION 



The random nature of ocean waves is most often descrited 
using the linear spectral model, and many engineering prob- 
lems related to offshore and coastal structures are solved 
by employing the concept of directional wave spectra. Eaves 
in the sea, however, can exhibit nonlinear properties. A 
departure of surface elevation from the Gaussian distribu- 
tion with positive skewness is an example. Phase velocities 
of high frequency component waves not satisfying the linear 
dispersion relation is another example. Breaking of waves 
is a spectacular example of nonlinear behaviour of water 
waves. 

The study of nonlinear, water waves traveling over a 
horizontal bottom dates back to Stokes ( 1 847) , who solved 
this problem in the form of a power series in terms of a 
small parameter ' £ * related to the average wave slope ( £ = 
ak, where a is the wave amplitude and k is the wavenumber) . 
It was found, however, that the Stokesian type series are 
nonuniformly convergent, and are only valid in deep and 
intermediate water depth. 

Boussinesg (1877) and Korteweg and de Vries (13S5) 
developed the cnoidal wave theories, which are based upon 
power series in terms of wave height relative to water 
depth. These power series are uniformly convergent in 
shallow water, but are invalid in deep water. 

Making the connection between the two above theories, 
Goda (1983) developed an empirical parameter to describe the 
phenomena of wave ncnlinearities with good results. The 
parameter bridges the wave steepness in deep water and the 
Ursell’s parameter Or = (a/d)/(d/L) 2 (with a the wave 
amplitude, d the water depth, and L the local wave length) 



in very shallow water, spanning the full range of water 
waves. 

When confronted with experimental results, Le Hehaute et 
al. (1968) demonstrated that the nonlinear theories were not 
much tetter than the linear wave theory; nonetheless there 
is a need to clarify the extent of wave nonlinearity in the 
sea so that engineering problems can be solved with much 
more confidence and accuracy. 

During the last two decades, numerical approaches based 
on truncated Fourier expansions have been proposed and well 
verified experimentally. The first such approach was devel- 
oped by Chappelear (1961) involving the use of the velocity 
potential. Dean (1965) used the stream function to develop 
his numerical wave tteory, which was computationally simpler 
than Chappelear’ s technique. Dean’s stream function form 
satisfies the Laplace equation, the kinematic free surface 
boundary condition, and the bottom boundary condition; the 
parameters in the stream function expression are chosen by a 
best fit to the dynamic free surface boundary condition. Von 
Schwinf and Reid (1972) developed characteristic solutions 
to finite amplitude waves using the velocity potential. 
Cokelet (1977) extended the method originally developed by 
Schwartz (1974) to allow a very accurate calculation of the 
characteristics of water waves. The procedure involves 
expressing the complex potential solution in a Fourier 
series and represents the Fourier coefficients in terms of a 
perturbation parameter. Rienecker and Fenton (1980) used a 
finite Fourier series to describe waves by solving a set of 
nonlinear equations using Newton's method. Recently, le 
Mehaute et al. (1984) developed a parameterized solution 
which is valid for all practical ranges of values of wave 
amplitude, frequency and water depth. In view of its rela- 
tive simplicity, and experimental validity, the parameter- 
ized solution is recommended for engineering applications. 
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The objective of this work is to utilize the spectral 
model by Le Mehaute et al. (1984) to compare with wave data 
acquired at Santa Barbara during the 2-6 February, 1980. 
The range of application of the model will be checked 
against measured wave spectral components and individual 
waves. The wave nonlinearity and Ursell number are analyzed 
using the model. Comments and conclusions are made for 
further study. 
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II. THEORETICAL BACKGROUND 



A. PARAMETERIZED SOLUTION TO THE NONLINEAR HATE PROBLEM 
1 . Ov ervie w 

The higher-order Stokian wave theories become alge- 
braically very complex. For practical purpose, it was 
desirable to develop wave theories that could be computed to 
any order. Le Mehaute et al. (1984) developed a parameter- 
ized solution to the monochromatic, irrotational , nonlinear 
water wave over a horizontal bottom. The formulation is 
expressed in a closed form in terms of a truncated Fourier 
series. All coefficients are expressed in terms of simple 
analytical functions which contain three parameters, namely, 
wave height, H, wave period, T, and water depth, d. The 
number of terms in the Fourier series is proposed paramet- 
rically for a given level of accuracy. The solution is 
assumed to be applicable for all practical ranges from deep 
to shallow water waves, i.e., d/Lo > 0.005. The deep water 
solution for d/Lo = 0. 5 is assumed to be valid for d/Lo > 
0.5. 



2. The oretic al Form ulat ion 

As in the periodic water wave problem the formula- 
tion has linear and nonlinear boundary conditions and a 
linear governing differential eguation. 

Choosing a coordinate system moving with wave speed 
c = (J / k so that the time dependency is removed, i.e., a 
stationary problem , it is possible to assign a constant 
value to the stream function at the free surface. The vari- 
able (J is the wave angular frguency ( (J =27 T /T) and k is 
the wavenumber (k = 2 7T / L) . The solution is a velocity 
potential function of the form 



0 { (x- ct),z} or 0{$,z) which satisfies the Laplace 
equation : 

y 2 0 = 0 (2.1) 

The variable Q = kx - (J t denotes phase, and the z axis is 
vertical positive upwards from the S.W.L. (still water 
level) and the horizontal x axis is positive in the direc- 
tion cx wave propagation. 

The kinematic bottom boundary condition is written 
as 





0 



( 2 . 2 ) 



Ike dynamic free surface boundary condition, is 
derived using Bernoulli’s equation written in the dimen- 
sionless form at z = Tj 



77 







\ 




( 2 . 3 ) 



where H is the wave height, JJ is the water wave surface 
elevation and Be = an arbitrary Bernoulli constant. 

The kinematic free surface boundary condition is 




( 2 . 4 ) 



The e's are the allowable errors which should be small. 

The measure cf how well these two boundary condi- 
tions are satisfied is defined by E j , and , which are 
the mean squared errors of the dynamic and kinematic free 
surface boundary conditions : 
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(2.5) 




( 2 . 6 ) 



in which j = sample at evenly spaced phase angles (e.g., 
with 5 ° interval) along the wave profile. The Bernoulli 
constant. Be, is determined such that the variation cf the 
phase speed c is minimized over a wavelength. For an exact 
solution, E j and E ^ wculd be zero. 

A solution in the form of a limited series (or trun- 
cated Fourier expansion) is assumed 




(2.7) 



in which (p = dimensionless potential function; An = a 
derived coefficient; n is an integer and N = the number of 
terms. All the internal wave field terms (velocity compo- 
nents, dynamic pressure and acceleration component’s) can be 
directly determined from this equation. Eg. 2.7 satisfies 
both the Laplace equation (2.1) and the kinematic bottom 
boundary condition (2.2) exactly. 

The equation for the free surface correpcnding to 
eg. 2.7 is: 



2 a 



K 

2 a 



N 



k X"* A n sinh I" nk ( d + 

2(J 2 L ^ 1 



n-l 




cos n 0 



( 2 . 8 ) 



in which K = an integration constant satisfying 




= 0 



(2.9) 
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At the wave crest ( Q = 0) , Tj - T} c and the wave trough 

( 0 = TT ) , T) = T} t -Therefore 



Ve + I’M = , + 

2 a 



( 2 . 10 ) 



in which E^ is to be determined and must be small. 

Ie Mehaute et al (1984) defines the phase velocity c 
or wavenumber k in accordance with the Hamiltonian varia- 
tional principle of minimum energy, which, in practice, is 
obtained through minimizing Zj with respect to the waven- 
umber, k 

d e, 

d * 0 i 2 -") 



Such a system of equations can eventually be solved as a 
function of phase velocity. 

The solution to the problem consists in providing an 
analytical solution which describes all the wave character- 
istics as functions cf only three parameters : wave height, 
H, wave period, T, and water depth, d. These parameters are 
grouped into two dimensionless parameters 



Depth parameter: 




d a ' 



2 7TS 



( 2 . 12 ) 



Wave steepness parameter: 



H_ _ o (jJ_ 
L* 7T 3 



(2. 13) 
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in which lo is the linear deep water wavelength and a = H/2 
is the wave amplitude. 

Having specified the initial boundary value problem 
and analytical forms cf egs. 2.3 and 2.4, the problem now 
consists of expressing all unknowns ; An, k, K, and N (eg. 
2.8) , in terms of the two parameters (5 and y , such that 
a desired level of accuracy is achieved. These unknowns are 
determined by satisfying the boundary criteria as closely as 
possible. 2g. 2.11, which specifies the phase velocity and 

wave number, must also be satisfied. The method of solution 
is one cf trial and error involving multiple iterative 
processes as expressed in the flow chart in Figure 1. 

For practical purposes, the final formulation may 
be presented in the following: 

INPUT DATA : Given E, T, and d then 

L o=gT 2 j2lT y **=27 T/z.* y k-2 7T I L ^ (J = 27T / T 

8 =*/u C 7 = h / l o p m y / t . lll , f4 ' 5 8 8 < 5 °'" 5 ; 

and (J 2 = $ k* tanh k*d . 



(The symbol * denotes the linear solution) . 
FORM CF SOLUTIONS: 




£ S 

a 




cosh 



\^n k ^d + j sin n 0 



(2. 14) 




a_S_k 

(J 2 



N 



I 

n- / 



Ap 



sink + C os nQ 



(2. 15) 
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INTEGRATION CONSTANT: For application purposes, 

ter ized form of K is determined by 



a parame- 



K = - /. 764 a_J_k 

(J 2 



o.55 o .77 

P 5 



(2. 16) 



The use of K as expressed above generally yields larger 
values cf E’ s than the value of K obtained by eg. 2.9 
directly. 

NONLINEAR WAVE NUMBER : k = 2 7T/L 
1. 0.005 < § <0.1 



k = k'= 



[ , + l.879p l -° 97 -0.094/]^ B j„ ( 2 X ,0^ 1 } 



0.528 



(2. 17) 



2. 0. 1 < (5 < 0.5 



k = k = 



/+ 3.466 fi 1 - 5 ft 0 ' 3 ' 5 1 



tanh k*d 



(2. 18) 



The solution valid for (5 = 0*5 is also considered valid for 
5^0.5 provided that all the parameterized functions which 
follow are determined with (5 = 0.5 irrespective of their 
real value; this means that deep water solutions corre- 
sponding to 5 > 0.5 are not influenced by water depth. 

COEFFICIENTS: 

Coefficient A j : 

A , = A*( / - R ) 

(2. 19) 
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( 2 . 20 ) 



cosh k*d 

1. 0.005 < 5 <0.1 

✓ 

* = «; (s-BSip)"' 



in which 



2 . 



r 3 r ^ 2,82 

a» y = 0.8132 exp (-0.0153 jin (1 0 J 0 )] } 

b ^ = 0.0747 exp (1. 1057 [in (ID 3 5 } j 75 } 
0.1 < 5 ^ 0.5 



= 0.3-i- (-0.3 + a" ) exp. 



bj[ln(l0 5)] Cl 



in which a" 

I 






c 



II 

/ 



17.30/J 2 ' 38 

-/ n /.836 

1.369 tanh '{24.09/J ) 

n 3.04 

1.3676 exp (- 1 79 . O/j ) 



Coefficient : 

1. 0.0C5 < 5 < 0. 1 



in 



in 



which 



= 0.497 A f exp[ a‘ 2 (5.882p ' 

a ^ = -0.00933 exp{1.494 [ln(10 3 8 )]°‘ 79 1 
b* 2 = -0.953 exp {-0.029 [ lri (10 3 5 )] 2 ' 0/ } 



2. 0. 1 < 5 < 0. 0 

// 

A 2 = 2.5x/ 0 5 / 4 / fxpjc* 5^] 

-0-95 

which a'^ = 10. 547 exp (-0.0134^ ) 



( 2 . 21 ) 



( 2 . 22 ) 



(2.23) 



(2.24) 
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_-/.7g 

b» 2 = -16 + 15. 318 exp (-0.0017/J ) 

_ / 

c n 2 = 5-4.328 exp (-0.0112^} ) 

Coefficient A j : 

i. o. J05 S 4 o. i 



<3 = 0.672 A 2 txp^a' 3 (5.2%2fl) b 3^ 



>.78 



in which a*^ = -0.021 exp{1.355 |ln(10 3 6 )j } 

/ 2 8 

b' 3 = -0.859 exp {-0 .0747 [ln(10 3 5 )] ' J 
2. 0. 1 < (5 4 : 0.5 



(2.25) 



exp ( c* ) 



^o.52 -o.4l8 

in which a” = -3.537 0 (j 



Coefficient ftp (n > 3) 



(2. 2c) 



r^o.52 _ 0 .4I8 

i /7 = **P (~ 3-537 0 ft ) 



(2.27) 



3. An a lyti c al Vali dity 

Ihe validity cf this theory is defined through the E 
values as derived by egs. 2.3, 2.4, 2.5, 2.6, 2. 1C and 2.11 
by order of importance. The magnitude of the E's depends on 
the number of terms N used in the series. 3 eyond a certain 
value of N, it is found that the E values tend towards 
asymptotic values, indicating that they are not reduced by 
adding more terms, i.e. , the series is bounded. With this 
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property in the solutions, the number of terms N reguired 
can be expressed in parameterized form by : 



N = integer part of {a^+ j 


J - 0.2)} 


(2. 28) 


n ; = 1 ♦ 19.13 exp {-0.344 £ln ( 1 0 3 l 


) )} c 1 - 12 >3 


(2. 29) 


n 2 = 32.72 exp{-0.212 jin (10* 5 )]< 


1 . 31 )} 


(2. 30) 


le Mehaute et al, (1984) showed that the parameter- 
ized solution approximates the free surface dynamic boundary 
condition as well as, if not better, than any other existing 



theory. It also approximates the kinematic condition gener- 
ally letter, except for the stream function theory devel- 
oped ly Dean (1974) by virtue of its definition. 

Iruncated Fourier series do not represent limit 
waves well since the cusp at the wave crest reguires an 
infinite number of terms as was demonstrated by Dean (1S74). 
Therefore, the accuracy of the parametric solution decreases 
rapidly as the wave heights approach limit wave conditions, 
which is where Ur is large. 

B. MEASURE OF N0H1IHEARITY 

The parameter which represents the wave nonlinearity in 
deep water is the wave steepness, or a measure of the 
surface elevation slcpe. The most common representation is 
either the ratio of wave height to wavelength H/L, or the 
parameter ak. 

In the region of shallow water waves, the Stokian wave 
profile for higher orders is an expansion using the nondi- 
mens icnal term 

Ur = (a/d) / { d 2 /L 2 ) (2.31) 
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as the perturbation parameter. This term is called the 
Orsell parameter and has been shown to govern the transfor- 
mation of water waves in shallow water. For example, it has 
been demonstrated that the nonlinear shoaling of water waves 
is predicted as the function of the Ursell's parameter. 

When Ur <<1/ the linear small amplitude wave theory 
applies. In principle, more and more terms of the power 
series are required in order to keep the same relative accu- 
racy as the Ursell parameter increases. 

The Ursell parameter is a usefull guide, but is not 
necessarily sufficient for judging the relative importance 
of the nonlinear effects. A qualitative idea of the impor- 
tance of nonlinearity is also given by the ratio of spectral 
values of each harmonic frequency to the spectral value of 
the fundamental frequency. However, the nonlinear effects of 
shoaling transformation are cumulative, so that the ncnli- 
nearities may be locally weak, but significant cross- 
spectral transfer can occur if the shoaling region is wide 
in comparison to the nonlinear interaction distance. 
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III. EXP ERI MENTAL PEOCEDOEES 



Save and velocity data acquired at Santa Barbara, 
Califcrria, during the 2nd, 3rd, 4th, 5th and 6th February 
1980, as part of the Nearshore Sediment Transport Study, are 
used to be compare with the parameterized model. The field 
data, considered in this work, include only measurements 
made by the current meters outside the breaking zone. 
Current meters C07X, C03X and C91X (see Figure 2) are used 
on the 2nd and 6th of February; current meters C03X, C01X 
and CODX are used on the 3rd of February, and CO IX and CODX 
are used on the 4th of February. 

For each day and each instrument the following statis- 
tics were calculated: average wave height, fT, the maximum 

wave height, Hmax, the significant wave height, Hl/3, the 
root mean squared wave height, Hrms, and the average of the 
heights cf the 1/10 highest waves, HI/10, as well as the 
associated water depths, and the surface 
elevaticn spectrum. 

A. HAVE SPECTRAL ANALYSIS 

The surface elevation spectra were calculated frcm 63 
minute current meter records. The water particle velocity 
components measured by the current meters were recorded on a 
special receiver/tape recorder, digitally low-pass filtered 
and sampled at 2 samples/sec. The data were then high-pass 
filtered at 0.05 hz to exclude low frequency variations. 
These current data were used to infer wave heights. The 
velocity signals were convolved using linear wave theory to 
obtain surface elevations. The complex Fourier spectra of 
the horizontal velocity components were first calculated and 
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vectorially added to give V (f) . The complex surf ace eleva- 
tion spectrun, X (f ) , was calculated applying the linear wave 
theory transfer function H(f) 

X (f ) = H(f) * v (f ) (3.1) 
where 

H (f ) = (sinh kd) / {(J cosh k(d + z)} (3.2) 

with z = measurement cf elevation. 

The calculated energy density spectra are to be compared 
with the parametric model. Since the parametric model 
defines all the truncated Fourier series coefficients as 
simple analytical functions in terms of the wave height H, 
wave period T, and water depth d, it was necessary to trans- 
form the energy density spectra of the sea-surface elevation 
into a discrete (spikes) spectra as shown in Figure 3. 

The discrete spectra- were calculated by first finding 
the fundamental freguency defined as the frequency associ- 
ated with the highest energy density peak. The harmonics 
were defined as the integer mutiples of this fundamental 
freguency. The bandwidths associated with the energy peaks 
were defined at the frequencies corresponding to the half 
power pcints of the energy density values. The variance of 
the peaks, calculated as the area of the energy density 
spectrum between half power points, describe the height of 
the spectral spike. The analysis was carried out to the 
highest harmonic for which half power points could be 
defined. The reason for defining the bandwidths at the half 
power points is because the half power values correspond to 
the level at which the spectral components are independent 
of one another. It was noted that the bandwidth increases 
approximately with the number of harmonic. 
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Assuming Tj - 0, the variance of the wave surface is 

given by 



var 





(3.3) 



The parameterized variance, or energy, can be defined from 

equation 2.15 in the discrete mode 

f 1 2 



var (n) = 







a 3 k 

a 2 






(3.4) 



with n = 1,2,3. ..N. 

The wave height (and wave amplitude) used to calculate 
the parametric energy line spectrum was the average wave 
height, H, because it was the wave height parameter which 
offered best fitting to the field data. 

Comparisons for all days were made between the converted 
actual line spectrum and the parametric line spectrum calcu- 
lated using the input parameters of average wave height H, 
average water depth d, and the period corresponding to the 
defined fundamental frequency (T = 1/f) , calculated for the 
63 minute measurements. 



E. INEIVIDUAL WAVE AKALYSIS 

A second approach was to compare individual waves meas- 
ured in the field with the parametric model predictions. 
The filtered and digitized velocity signals were convolved 
using linear theory tc obtain surface elevations. The indi- 
vidual waves were defined using the zero-up crossing method 
as explained in Figure 10. The highest maximum and the 
lowest minimum of the surface elevation within a period 
interval define the crest and the trough of a wave. A wave 
height Hi is defined as the range of TJ (t) in that 
interval, the period is the time interval between two 
consecutive zero-up crossings cf T) ( t) . 
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A Fourier analysis was performed on each individual wave 
using a Fast Fourier Transform algorithm. The parametric 
coefficients associated with the surface elevation formula 
were also determined for each wave using the measured Hi, 
Ti, and di. Statistics of Fourier coefficients from field 
data and the parametric model were generated for each day 

N u 

Tn = (1/Nw) * ^An(i) (3.5) 

/=/ 

(where Nw = number of analyzed waves) . For each run compar- 
isons were made between this approach and the wave spectral 
analysis. Since the average wave period was about 15 
seconds, the total number of waves in each 68 minutes record 
was about 300, which give reasonable wave statistics. 

C. SUMMARY 

Beth analyses were based on the similarity between the 
Fourier series representation of the sea surface elevation 

CO 

77 (t) = aO + ^ an cos n Q (3.6) 

n=i 

(with aO = mean) and the sea surface elevation described by 
the parametric model as defined by the equation 2.15. The 
coefficients to be compared are the an’s for the Fourier 
series defined by the Fourier integrals and the Le Mebaute 
coefficients 

- (1/ CT **2) *a*g*k* An * sinh lnk[d.+ Tj)} 
where An are defined parametrically as function of the depth 
parameter (5 and wave steepness ~y ,and Tj defined itera- 
tively from equation 2.15. 
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The spectral data analysis is a linear approach. The 
spectral components are assumed independent of each other. 
The analysis considers that the spectral component phase 
angles at each freguency are random (putting the problem in 
a probabilistic or stochastic setting). On the other hand, 
the individual wave analysis is also a linear approach, but 
the phase relationship between harmonics is presumed deter- 
ministic (i.e. phase coupled with the fundamental freguency, 
as evidenced by the observed peaked crests and flattened 
troughs) . 
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IV. RESULTS 



A. SPECTRAL ANALYSIS 

1 . No nlin ear Statistics 

A typical energy density spectrum for all analyzed 
instruments is presented in Figure 2 (upper panel) . The 
presence of strong harmonics in the spectrum indicates the 
importance of nonlinearities of the waves in this region of 
shallow water. A qualitative idea of the importance of 
nonlinearity is given by the ratio of the harmonic energies 
to the primary frequency energy E{nf)/E(f), where n = 2, 3 
... are the harmonic frequencies. Due to the difficulty of 
defining the half power points for the harmonics higher than 
the third, only the ratios involving the first and second 
harmonics are presented. 

Tables I - V give the ratios for all five days 
studied for both the parame terized model (denoted L.M.) and 
the field data (denoted by F.).The water depth and the 
Ursell parameter associated with each measurement are also 
presented. The range of water depths is from about 670 cm 
to 138 cm. The range cf the Ursell parameter is from about 6 
to 12C. 

For large Ursell numbers ( Ur > 75 ) , the parameter- 
ized model overestimates the importance of the energy asso- 
ciated with the harmonics relatively to the fundamental 
frequency energy. The same ratios for the field data are 20 
- 60 percent lower. For small Ursell numbers ( Ur < 25 ) , 
the parameterized model underestimates the importance of the 
harmonics relatively to the fundamental frequency energy. 
The same ratios for the field data are 0.5 - 40 percent 
higher. 
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The best fitting between the parameterized model and 
the field data occurs at intermediate Ursell numbers ( 25 < 
Dr < 75 ) . It is found that the parameterized model energy 
ratios compared with the field data energy ratios have an 
error of less than 20 percent for the cases of intermediate 
Drsell number. Some exceptions were observed to the above 
indicated percentages, but the energy contribution at these 
frequencies is very small. 

2. local S pectral Comparisons 

An absolute comparison was made between the parame- 
terized model spectral energy peaks and field data spectral 
energy peaks as a function of frequency. Typical results for 
three different Drsell numbers are presented in Figures 4 
and 5. In regions of strong nonlinear effects ( large 
Drsell number) , the parameterized model underestimates the 
energy at the fundamental frequency. The energy drops 
linearly from the fundamental frequency to the harmonics, in 
the parameterized model, while field data spectra fall off 
much faster. In regions of small Drsell numbers the parame- 
terized model concentrates almost all the energy in the 
fundamental frquency, overestimating it relative to the 
field data. The best agreement between the parameterized 
spectral energy peaks and field data energy peaks are 
observed for intermediate Drsell numbers. 

The Fourier coefficients used to evaluate field data 
spectra are plotted against the Le flehaute spectral coeffi- 
cients ( 1 / (J **2) agkA sinh[nk(d +7^)} (Figs. 6-8). A 

reasonable correlation is found for the first three 
coefficients. 

The spectral analysis approach is summarized by 
plotting the ratio of the Le Mehaute spectral peaks to field 
data spectral peaks against the Drsell parameter (Fig. 9) . 
The ratio for the fundamental frequency energy decreases as 
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the Orsell parameter increases emphasizing how the parame- 
terized model underestimates the fundamental frequency spec- 
tral peak in situations of strong nonlinearities ( large 
Ursell parameter) . The ratios for both the first and second 
harmonic frequency energies increase as the Orsell parameter 
increases; it must be noted, as indicated earlier, that the 
data only ranges in Orsell numbers from about 6 to 120. 

B. IHBI7ID0AL WAVE 2HA1YSIS 

Individual waves were identified using the zero-up-cross 
method. A Fourier analysis was performed on each individual 
wave. Ihe parameterized spectral coefficients were also 
determined for each individual wave. 

The Fourier coefficients associated with the fundamental 
frequency, first and second harmonics are plotted separately 
against the corresponding parameterized coefficients. The 
coefficients at the fundamental frequency (Fig. 14) are 
reasonably well correlated (correlation coefficient = 
0.80). The Fourier coefficients plotted against the Le 
Mehaute coefficients follow a straight line with intercept 
at the origin and slcpe of abcut one. The group of plotted 
points with a greater slope (approximately 2) are associated 
with measurements in shallower waters with higher Ursell 
numbers (see Figs. 12-16). The graphs of the the first ar.d 
second harmonic coefficients (Figs. 17-18) show similar 
char acteristics as the fundamental frequency; however, the 
correlation between the coefficients is poorer. 

The ratio of the le Mehaute coefficients to the Fourier 
coefficients is also plotted against the Ursell parameter 
The ratio for the fundamental frequency coefficients 
converge to a bound value of 0.5 at large Ursell numbers 
(Fig. 19). For Ursell numbers less than 75 the ratio of the 
le Mehaute to Fourier coefficients is about 1 with the 
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exception of a few points near Dr = 0 which were identified 
to be related with noisy signals. For large Orsell numbers 
the Fourier coefficients are about twice the Le Mehaute 
coefficients, again showing the tendency that the parameter- 
ized model has to underestimate the fundamental peak energy 
in regions of strong nonlinearities (large Ursell numbers). 

The ratios for the first and second harmonic coeffi- 
cients (Figs. 19 and 20) show a bounded value near unity as 
the ursell parameter increases. Noisy signals near Ur = 0 
are also observed. 

C. COMPARISONS BETWEEN APPROACHES 

Beth linear (spectral analysis) and nonlinear (indi- 
vidual wave analysis) approaches lead to similar results. 
The parameterized model is observed to overestimate the 
fundamental frequency energy in regions of small Orsell 
number (weak nonlinearities) and underestimate it in regions 
of large Ursell number. The individual wave analysis pres- 
ents a limit tendency for large Ursell numbers, with the 
fundamental frequency Fourier coefficients about twice those 
of the parameterized spectral coefficients. The better 
correlation observed between the Fourier coefficients and 
the le Mehaute coefficients in the individual wave analysis 
as compared with spectral analysis approach is due, perhaps, 
to the definition of half-power criteria or to the averaging 
processes involved in calculating the wave spectra. 

The reason for the similarity of results for the two 
approaches can be seen by comparing the spectra calculated 
in the usual manner with a "nonlinear spectrum". A nonli- 
near spectrum was calculated using the Fourier amplitude 
coefficients calculated for the individual waves. Using the 
same bandwidths as the ordinary spectrum, amplitude coeffi- 
cients were sorted into the respective frequency bands and a 
nonlinear spectrum calculated as 
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Gnl (f) 



/ 24f 



(4.1) 



= (1/Nw)^A?(f) 

i=! 

where the coefficients Ai(f) are centered on frequency f 
within bandwidth A f. The total number of waves is denoted 
by Nw. An example of the comparison of the nonlinear spec- 
trum calculated in this manner with the linear spectral 
method is shown in Figure 21. The two spectra are similar 
and bcth show strong harmonic peaks. The spectral energy in 
the ordinary spectrum tends to be smeared more, whereas the 
peaks and valleys in the nonlinear spectrum are more sharply 
defined. 
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?. CONCLUSIONS 



The coefficients of the parameterized solution tc the 
monochromatic, irrotational, nonlinear wave over a hori- 
zontal bottom, as developed by Le Mehaute et al. (1984) , are 
expressed in terms of simple analytical functions of wave 
height, H, wave period, T,and water depth, a. All quanti- 
ties explicitly derived from the potential function 
(velocity components, dynamic pressure, acceleration, i.e. , 
the internal wave field) can be directly determined. The 
relative simplicity of the calculations of the parameterized 
solution recommends it for engineering applications. 

The wave spectra determined from the parameterized solu- 
tion show good agreement with the measured wave spectra from 
field data for conditions of intermediate Ursell numbers (25 
< Ur < 95) . For small Ursell numbers, the parameterized 
solution concentrates almost all the energy in the funda- 
mental frequency representing a sinusoidal waveform. Under 
these conditions, the field wave spectra always contained 
relatively less energy associated with the fundamental 
frequency and more energy distributed to the harmonics. 

For large Ursell numbers, the parameterized solution 
describes a waveform with more peaked crests and flatter 
throughs relative to the measured waves, since it underesti- 
mates the energy associated with the fundamental frequency 
and overestimates the energy associated with the higher 
harmonics. 

The ratios between the Le Mehaute and field data spec- 
tral coefficients associated with the 1st and 2nd harmonics 
present a bound limit tendency of 1 for large Ursell 
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numbers, indicating good agreement between the two processes 
for these conditions. The same ratios for the coefficients 
associated with the the fundamental frequency show a bound 
limit tendency of 0.5 , indicating that the first coeffi- 
cient of the parameterized model needs to be corrected for 
strong (Or > 75) nonlinear effects. 

Further comparisons should be made for different sea 
conditions to determine the corrections to the parameterized 
model coef f icients . However, it must be noted that a limi- 
tation of the parameterized solution is that its accuracy 
decreases rapidly as the wave heights approach limiting wave 
conditions due to the use of a truncated Fourier series. 
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Figure 1 Flow Chart of the Parameterized Nonlinear 

Wave Model. 
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SANTA BARBARA BEACH 

CALIFORNIA 

2-6 FEB. 1980 



WAVE STAFFS 




Figure 2 Beach Profile with the Instruments Location. 
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frequency (Hz) 




Figure 3 Typical Energy Density Spectrum (above) and 
Conversion into Discrete Energy Spectrum (below). 
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TABLE I 



Batio of Harmonic Energy to Primary Frequency Energy 

02 FEB- 80 



DEPTH (cm) 


212.5 


333.7 


398.7 


669.6 


URSELL No. 


115.6 


31.3 


23.4 


7.6 




L.M. 


E. 


L.M. 


F. 


L.M. 


F. 


L.M. 


F. 


E(2f)/E(f) 


0.548 


0.165 


0.167 


0.132 


0.108 


0.120 


0.008 


0.021 


E(3f )/E(f ) 


0.236 


0.056 


0.019 


0.061 


0.008 


0.028 


- 


- 


Instrument 


C07X 


C03X 


C01X 


C0DX 



TABLE II 

Batio of Harmonic Energy to Primary Frequency Energy 

03 FE E. 80 



DEPTH (cm) 


302.3 


379.8 


650.7 


UBSELL No. 


53.4 


31.7 


11.1 





L.M. 


F. 


L.M. 


F. 


L.M. 


F. 


E(2f )/E(f ) 


0.365 


0.194 


0.173 


0.150 


0.025 


0.145 


E(3f)/E(f) 


0.097 


0.073 


0.021 


0.064 


0.008 


- 


Instrument 


C03X 


C01X 


C0DX 
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TABLE III 



Ratio of Harmonic Energy^to Primary Fre< 3 uenc Y Energy 



DEPTH (cm) 


378.8 


649.7 


URSELL No. 


31.2 


7.9 


_ 


L.M. 


F. 


L.M. 


F. 


E(2f)/E(f) 


0.171 


0.159 


0.013 


0.125 


E(3f )/E(f ) 


0.020 


0.054 


0.001 


0.014 


Instrument 


C01X 


C0DX 



TABLE I¥ 



Eatio of Harmonic 



:rgy to Primary 
05 FEE- 80 



Frequency Energy 



DEPTH (cm) 


363.6 


634.6 


URSELL No. 


24.8 


6.8 




L.M. 


F. 


L.M. 


F. 


E(2f )/E(f ) 


0.123 


0.101 


0.008 


0.069 


E(3f)/E(f ) 


0.01 


0.050 


0.001 


0.034 


Instrument 


C01X 


C0DX 
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TABLE ? 




DEPTH (cm) 


138.8 


278.7 


345.4 


URSELL No. 


115.5 


23.8 


11.5 


— 


L.M. 


F. 


L.M. 


F. 


L.M. 


F. 


E(2f )/E(f ) 


0.548 


0.186 


0.113 


0.196 


0.027 


0.110 


E(3f )/E(f ) 


0.236 


0.068 


0.009 


0.088 


0.002 


0.049 


Instrument 


C07X 


C03X 


C01X 



LOCAL SPECI HAL COMPARISONS LOCAL SPECTRAL COMPARISONS 




L> 

I 

o 

03 

23 

U 1 

L*, 

CV 

<D 



o 

2: 

-j 

-3 

id 



r 3 



0 021 0 011 0 001 



0 09 0 04 0 09 0 0? OOt 

(2..KD) - U)ADH3N3 




o n 0 01 oot 0 00 00* 



0 0* O K 0 OC 0 fit 

(r-.HO) - CJJADH3N3 



42 



URSELL NO. = 11.5 6FEU80-C01X 



ENERGY(F) - (CM**2) 

0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0 110.0 120.0 



LOCAL SPECTRAL COMPARISONS 



0.00 



0.05 



o 

□ 



LEGEND 
□ = FFTCC 
o = LE ME1IAUTE 



0.10 0.15 

FREQ. (HZ) 



□ 

o 

T 

0.20 



0.25 



URSELL NO. = 31.3 2FEB80-C03X 



Figure 5 local Spectral Comparisons for 
Intermediate Orsell Number. 
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L.M. - F. COEFFICIENT RELATIONS 
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L.M. - F. COEFFICIENT RELATIONS 
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LM. - F. COEFFICIENT RELATIONS 
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SPECTRAL PEAKS COMPARISONS 
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where 



X is a crest 

* is a trough 

0 is the zero-up crossing 



Figure 10 Defining Haves using 
Zero-up Cross Method. 
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Figure 11 Coefficient Ratio at Fundamental 
Freguency against the Orsell Number. 
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Individual Waves 
2 Feb. 80 
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Figure 17 Correlation of the Fundamental Frequency Coefficients 

for the Current Meter C01 - 2Feb. 80. 
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Figure 18 Correlation of the Fundamental Frequency Coefficients 

for the Current Meter COD - 2Fefc. 80. 
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Individual Wave 
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Figure 21 Comparison between Normal Spectrum 
and Nonlinear Reconstructed Spectrum. 
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